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For electronic states and photoinduced charge dynamics near the neutral-ionic transition in 
the mixed-stack charge-transfer complex tetrathiafulvalene-p-chloranil (TTF-CA), we review 
the effects of Peierls coupling to lattice phonons modulating transfer integrals and Holstein 
couplings to molecular vibrations modulating site energies. The former stabilizes the ionic phase 
and reduces discontinuities in the phase transition, while the latter stabilizes the neutral phase 
and enhances the discontinuities. To reproduce the experimentally observed ionicity, optical 
conductivity and photoinduced charge dynamics, both couplings are quantitatively important. 
In particular, strong Holstein couplings to form the highly-stabilized neutral phase are necessary 
for the ionic phase to be a Mott insulator with large ionicity. A comparison with the observed 
photoinduced charge dynamics indicates the presence of strings of lattice dimerization in the 
neutral phase above the transition temperature. 

KEYWORDS; photoinduced phase transition, neutral-ionic transition, electron-lattice interaction, electron- 
molecular-vibration coupling 



1. Introduction 

Nonequilibrium electronic states have drawn attention 
as a field where novel electronic functions and properties 
are being searched for. Correlated electron and electron- 
phonon systems with rich electronic phase diagrams are 
promising because synergy can be exploited, which is in- 
herent to respective orders at low temperatures. Among 
nonequilibrium phenomena, photoinduced phase transi- 
tions occur synergistically since a low density of photons 
relative to the density of molecules or atoms can alter 
the electronic phase realized in such a system. A variety 
of photoinduced phase transitions are known to proceed 
very rapidly.^' 

From the viewpoint of controlling functions in nonequi- 
librium states, it is important to learn a means by which 
a photoinduced phase transition proceeds. The most fun- 
damental information is on relevant interactions that sta- 
bilize the ground state. It is generally true that, even if 
a few interaction parameters are sufficient to describe 
equilibrium properties, they are still insufficient for de- 
scribing nonequilibrium dynamics. Additional interac- 
tions can be crucial in determining whether and how the 
phase transition is photoinduced. Among organic mate- 
rials that have been regarded as strongly correlated elec- 
tron systems, some have been realized to have substan- 
tial electron-molecular-vibration (EMV) couplings that 
contribute to the stabilization of the ground state and 
photoinduced dynamics. 

For instance, the charge order in the quasi-two- 
dimensional organic salt a-(BEDT-TTF)2l3 [BEDT- 
TTF=bis(ethylenedithio)tetrathiafulvalene] is basically 
driven by a long-range Coulomb interaction."^' "'^ Rela- 
tively weak Peierls couplings to lattice phonons or molec- 
ular rotations are important for understanding the inter- 
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relation between the altered charge distribution and the 
structural deformation^' as well as the long-time behav- 
ior of photoinduced dynamics.^' From the early-stage 
dynamics after photoexcitation, however, substantially 
strong Holstein couplings to molecular vibrations are ev- 
ident in the Fano destructive interference between vi- 
brations containing C=C stretching and correlated elec- 
trons' motion.*^' These molecular vibrations stabilize the 
charge order to a substantial extent. 

Most organic compounds that show electronic phase 
transitions have C=C bonds in their constituent 
molecules. Their vibration frequencies arc well known to 
depend on the charged state of the molecule. It is there- 
fore expected that EMV couplings arc generally strong 
even though the electronic state has strong electron cor- 
relations in terms of intermolecular spin and charge cor- 
relation functions. In this study, we focus on the mixed- 
stack charge-transfer complex TTF-CA, whose thermal 
and photoinduced neutral-ionic transitions have inten- 
sively been studied both experimentally^°"^^^ and theo- 
retically.25-39) We review the origin of the discontinuities 
in ionicity and optical conductivity during the neutral- 
ionic transition. Their quantitative understanding will 
become important for controlling photoinduced phase 
transition dynamics. 

2. Discontinuities during Neutral-Ionic Transi- 
tion 

The importance of the intcrsite Coulomb repulsion V 
for discontinuities during the neutral-ionic transition has 
been pointed out in ref. 25. Without coupling to lattice 
phonons, i.e., without dimerization, the spin gap at the 
phase boundary is zero on the ionic (i.e., regular-Mott- 
insulator) side and V on the neutral (i.e., band- insulator) 
side in the limit of vanishing transfer integrals. Here, 
the creation of a neutral-ionic domain wall requires V/2 
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of energy, which ahows for the phase separation of neu- 
tral and ionic domains at the boundary. The transfer 
integral to introduces quantum fluctuations. As long as 
2ta is smaller than V, however, the creation energy re- 
mains finite. The discontinuous neutral-ionic transition is 
demonstrated by quantum Monte Carlo simulations. 

The ionic phase is known to have dimcrization. The 
strongest candidate for its mechanism is due to the spin- 
Peierls instability. A regular ionic phase is a one- 
dimensional half-filled paramagnetic Mott insulator, so 
that it can be described by a spin-1/2 antifcrromag- 
netic Hcisenbcrg chain. If the supcrcxchangc interaction 
is modulated by lattice phonons, the system is stabi- 
lized by dimcrization that produces a finite spin gap. The 
dimcrization in the ionic phase is demonstrated by quan- 
tum Monte Carlo simulations.^^) As lattice phonons are 
strongly dimerized, i.e., as the electron-lattice coupling is 
strengthened, the phase boundary is shifted because only 
the ionic phase is stabilized, which decreases the ionicity 
on the ionic side of the phase boundary. The discontinu- 
ity in ionicity is strongly reduced by the electron-lattice 
coupling. As a consequence, the character of the ionic 
phase as a Mott insulator is weakened, and the transi- 
tion becomes close to a continuous Pcicrls transition. 

For the discontinuous ionicity change, Painelli and 
Girlando have pointed out that both EMV couplings and 
intersite Coulomb repulsion V are important and that 
they collaborate. For finite systems, they perform ex- 
act diagonalization in the limit of a large on-site repulsion 
U and the second-order perturbation theory with respect 
to the electron-lattice (i.e., Peierls) and EMV (i.e., Hol- 
stcin) couplings. The electron number configuration is 
roughly described as 2020 in the neutral phase and as 
1111 in the ionic phase, so that the Holstein couplings 
stabilize the neutral phase, which are in contrast to the 
Peierls coupling that stabilizes the ionic phase through 
dimerization. 

As noted above, both of the Peierls and Holstein cou- 
plings are suggested to be necessary for the reproduction 
of ionicity and optical conductivity near the neutral- ionic 
phase boundary. In spite of this, the Holstein couplings 
have often been ignored for the following reason. Without 
the Peierls coupling, the bond-charge densities cannot 
be dimerized, and the ionic phase cannot have the fer- 
roelectric order. In contrast, the Holstein couplings only 
enlarge the site energy difference between the neighbor- 
ing orbitals in the neutral phase, so that they can be 
absorbed into the renormalized site energy difference in 
the adiabatic limit. The situation is quite similar to a 
phase transition described within the mean field theory 
in the sense that the Holstein couplings stabilize the neu- 
tral phase in a synergistic manner and enhance the dis- 
continuity during the transition. Namely, the degrees of 
renormalization are different between the two phases. 

In this paper, we include both couplings in the model. 
The charge dynamics after photoexcitation of the neutral 
phase shows a large contribution from molecular vibra- 
tions. This is consistent with the ultrafast charge and 
molecular-vibration dynamics recently observed by the 
transient reflectivity measurements. ■^'^^ Numerical calcu- 
lations show that the initial neutral state must be ac- 



companied with a small dimerization to reproduce the 
experimentally observed, coherently oscillating dimcriza- 
tion immediately after photoexcitation. The presence of 
strings of lattice dimcrization in the neutral phase has 

indeed been observed by the latest X-ray diffuse scatter- 
ing. 23) 

3. Model with Lattice and Molecular Phonons 

We use the one-dimensional half-filled extended ionic 
Hubbard-Peierls-Holstein model 
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creates an electron with spin a at site j, 

and nj='Y^^nja. The parameter denotes 
the transfer integral on a regular lattice (i.e., without 
lattice distortion), A the site energy difference between 
neighboring orbitals when molecular distortions are ab- 
sent, U the on-site repulsion strength, and V the nearest- 
neighbor repulsion strength. The lattice displacement Uj 
at site j modulates the transfer integral between the 
(j — l)th and jih orbitals and that between the jth and 
{j + l)th orbitals with the coefficient =Fa. The displace- 
ment wj"*^ in the mth mode on the jth molecule modu- 
lates the site energy with the coefficient fi'f^^ . The quan- 
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tic coefficients, and oJq. and uj^^ 
energies, respectively. 

For these model parameters, we take eV as the unit 
of energy and use <o=0.17,'^'') U^l.b, and F=0.6; we 
vary A around the phase boundary depending on the 
Peierls and Holstein coupling strengths. We define the 
strengths of these couplings as Aq, = a^/Ka and 

\f}j = J2m l^j"^^^ / ^pj^ ■ ^ong as the displacements 
are treated classically, the ground state is given by their 
static configuration, and the ground state is determined 



not by the distribution of /?, 



(m)2 



but by their sum, 

(m)_ 



A^j.'*'^) The displacements are scaled using a=Pj'=l, 

so that we have Aq = 1/Ka and Xpj = 1/K^^\ For 
simplicity, we set \/32i-i = ^I32i = A^. 

As for phonons, we take one mode for the donor 
molecule and two modes for the acceptor molecule in 
addition to the lattice phonon mode, and use param- 
eters that approximately reproduce the experimentally 
observed phonon energies: a;Q=0.013 (this energy is soft- 
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cncd to be about 0.0066 when strongly coupled to elec- 
trons), Lu^^2i = '^^iAi=0.040, w^al-i = w^D=0.055, and 

(2) 

iOp2i = <^/3A2=0.12. Donor and acceptor molecules are 
specified by odd and even j's, respectively. For simplic- 
ity, we set -?^^2i~^^2i because the conclusion docs not 
depend on the detailed distribution of EMV-coupling 
strengths. The displacements uai, wd, and va2 corre- 
spond to CA's QgV^ (320 cm-i), TTF's flgi/g (438 cm"!), 
and CA's agV^ (957 cm"-'^) modes, respectively, in rcf. 24, 
whereas we ignore TTF's agVc, (740 cm~^) mode because 
its coupling to electrons is very weak. 

Photoexcitation is introduced through the Peierls 
phase 
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The time-dependent vector potential A{t) for a pulse of 
an oscillating electric field is given by 
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where Wpmp is the excitation energy, Tpmp is the pulse 
width, and F is the electric field amplitude. 

The time-dependent Schrodinger equation for the ex- 
act many-electron wave function on the chain of N^12 
sites with periodic boundary condition is numerically 
solved by expanding the exponential evolution operator 
with a time slice dt=Q.Q2 eV^^ to the 15th order and 
by checking the conservation of the norm.^*^) The initial 
state is set in the electronic ground state. The classi- 
cal equations for the lattice and molecular displacements 
are solved by the leapfrog method, where the forces are 
derived from the Hellmann-Feynman theorem: 
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Unless otherwise stated, the initial displacements are 
at the minimum adiabatic potential for the electronic 
ground state, so that the initial velocities are zero. 

4. Ground-State Properties 

Figures 1 and 2 show the ionicity 
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as a function of the site energy difference A near the 
phase boundary with different combinations of and 
A^. In the vicinity of the phase boundary, one phase is 
stable and the other is metastable. By comparing their 
energies, only the ionicity in the stable phase is plot- 
ted. It is clearly shown that, as A^ increases, the neu- 
tral phase is stabilized, the discontinuity in ionicity is 
enlarged, the ionicity in the neutral phase on the large- 
A side of the phase boundary becomes smaller, and the 
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Fig. 1. Ionicity p as a function of site energy difference A for (a) 
weak Peierls coupling Aa=0.05, and (b) strong Peierls coupling 
=0.167 with different strengths of A^. 



ionicity in the ionic phase on the small- A side becomes 
larger (Fig. 1). As a consequence, in order for the ionic 
phase to be a typical Mott insulator, i.e., in order for 
the ionic phase to have nearly one electron per site, Xp 
should be so large that the neutral phase is sufficiently 
stabilized. In particular, when Aq is large, A^ needs to 
be large [Fig. 1(b)]. Otherwise, the ionicity is too small 
even in the ionic phase. 

Figure 2 has essentially the same information as Fig. 1. 
It is evident that, as Aq increases, the ionic phase is sta- 
bilized, and the discontinuity at the transition is sup- 
pressed. When A;3 is large, the phase boundary does not 
shift monotonically as a function of Aq [Fig. 2(b) and 
A/3=0.33 (not shown)]. This is caused by the competi- 
— ' + c - ) 



tion between the bond charge X^o-i'-Jcr'^i+io- 

localized by Aq and the site charge ^^{(^ja^^jiy) localized 
by A^. The coupling Aq stabilizes only the ionic phase, 
while A^ stabilizes both phases (i.e., the neutral phase 
strongly and the ionic phase weakly). When Aq is small 
and the dimerization is weak (i.e., when lattice phonons 
only slightly lower the energy of the ionic phase), the in- 
creasing Aq, gains energy by localizing the bond charge, 
but it loses more energy by delocalizing the site charge. 
As a consequence, for large A^, the intermediate Aq sta- 
bilizes the neutral phase more strongly than the ionic 
phase. This counterintuitive result cannot be obtained 
using the second-order perturbation theory. 2'^) Neverthe- 
less, the shift of the phase boundary is small in this case. 

Experimentally, it is known that the thermal transi- 
tion from the neutral to ionic phases is accompanied by 
a large ionicity jump from 0.3 to 0.7. This indicates a 
large A^g. Of course, a finite Aq is certain to induce dimer- 
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Fig. 2. lonicity p as a function of site energy difference A for 
(a) weak Holstcin coupling A^=0.10, and (b) strong Holstein 
coupling A^=0.25 with different strengths of A^. 



ization and the three-dimensional ferroelectric order with 
a broken inversion symmetry. ^'^^ As shown below, a large 
Aq is also necessary for the reproduction of the optical 
conductivity alteration during the transition [Fig. 3(b)]. 
Therefore, both Ac and A^ are large in TTF-CA. 

The optical conductivity a{uj) in the ground state | tpo) 
is calculated using 
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where j is the current operator j = —dH/dA, e is a peak- 
broadening parameter set at 0.05, and Eq = (i/'o I I V'o) • 
Figure 3 shows ct{uj) in both phases close to the phase 
boundary. When only A^ is large, cr(w) is largely altered 
during the transition [Fig. 3(a)]. The optical gap in the 
neutral phase strongly stabilized by A^ is large, while 
that in the ionic phase is small in the present choice of U 
for TTF-CA. Experimentally, such a large alteration is 
not observed at the transition. When only Aq is large, the 
discontinuities in physical quantities at the transition are 
suppressed: the (j{uj) spectra in both phases are similar 
[Fig. 3(c)]. 

5. Photoinduced Dynamics from Ground State 

In the case of a large Aq and a large Xp [Fig. 3(b)], the 
neutral phase near the phase boundary is photoexcited 
with an energy a;pmp=0.65 just above the optical gap. 
The time evolution of the ionicity p{t) during and after 
photoexcitation is plotted in Fig. 4(a). For comparison, 
the displacement on the acceptor molecule VAiit) (with 
a lower oj^a) is shown in Fig. 4(b), the displacement on 
the donor molecule —VD{t) [(—1) is multiplied so as to 
oscillate in the same phase with p{t).] in Fig. 4(c), and 
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Fig. 3. Optical conductivity spectra in ionic (with smaller A) and 
neutral (with larger A) phases near the phase boundary with 
(a) Aq=0.05 and A^=0.33, (b) Aa=0.167 and A^=0.20, and (c) 
Aa=0.20 and A^=0.05. 



the displacement on the acceptor molecule vpj2{t) (with 
a higher w^a) in Fig. 4(d). As p{t) increases, the electron 
density increases for the acceptor molecule and decreases 
for the donor molecule [eq. (6)], so that the displacement 
increases at the acceptor molecule and decreases at the 
donor molecule [eq. (5)]. Thus, the quantities shown in 
Figs. 4(b)-4(d) basically behave as (—1) times the co- 
sine function. The ionicity pit) receives a positive feed- 
back from these molecular displacements and oscillates 
in the same phase with them. This —cosine behavior is 
consistent with the experimental observation,^**^ demon- 
strating that the neutral phase is stabilized by EMV cou- 
plings. 

To investigate quantum effects, we treated the molec- 
ular vibration with the highest phonon energy quantum- 
mechanically as in ref. 9. The difference between the wave 
profile when classically treated and that when quantum- 
mechanically treated was small (not shown). This is in 
contrast to the case reported in ref. 9, where the pho- 
toinduced transition is from an insulator to a metal, so 
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that low-energy electronic excitations exist and interfere 
quantum-mechanically with phonons. In such a case, the 
difference is large. In the present case, where the pho- 
toinduced transition is from an insulator to another in- 
sulator, the electronic excitations maintain a gap much 
larger than the phonon energies. The wavelet analysis 
does not show a trace of quantum interference. 

In the present calculation, we did not adjust the dis- 
tribution of /3]™^^/i<r^™^ so as to fit to the experimental 
data. Therefore, we can only roughly compare the cx- 
perimental^^^ and theoretical results for contribution of 
molecular vibrations to the photoinduced ionicity mod- 
ulation. They are comparable in that the oscillation am- 
plitude in the photoinduced ionicity change amounts to 
approximately one third in the present calculation and 
one fifth in the experimental result. We have calculated 
the photoinduced dynamics for different combinations of 
(Aa, A/3), (0.167, 0.2), (0.167, 0.25), (0.167, 0.33), and 
(0.2, 0.33), all of which roughly reproduce the disconti- 
nuities in ionicity and optical conductivity at the transi- 
tion. In these cases, the above contribution of molecular 
vibrations is roughly one third. 

However, there is a qualitative difference between the 
experimentally observed and present numerical results, 
which originates from the fact that the initial state here 
is the ground state with the regular and static lattice 
configuration. It does not contain thermal fluctuations. 
In the experiment, the lattice phonons are ready to be 
dimerizcd.^'*^ The photoinduced ionicity dynamics has a 
large contribution from their slow oscillation, which fits 
to a damped oscillator. In the present case, the inversion 
symmetry of the neutral phase is broken by the oscillat- 
ing electric field, so that lattice phonons begin to dimer- 
ize. However, its growth rate is quite low because it essen- 
tially corresponds to the spontaneous broken symmetry. 
The dimerization becomes significant only after t=600 in 
Fig. 4. If the initial state contains random numbers in Uj, 

Vj"^\ iij, and Vj"^\ according to the Boltzmann distribu- 
tion at a finite temperature of 0.01 cV, these fiuctuations 
accelerate the growth of {—lyuj, but they do not oscil- 
late it like (—1) times the cosine function as experimen- 
tally observed. It should be noted that this calculation 
cannot simulate thermally induced domains larger than 
the present system. 

6. Dynamics with Precursory Dimerization 

In this section, we introduce in the initial state a 
dimerization (— 1)^m-,=0.01, which is much smaller than 
(-l)-'uj=0.058 of the ground state at A=0.218 on the 
ionic side of the phase boundary. After obtaining the 
electronic ground state with fixed (— l)-'uj=0.01, we take 
it as the initial state and apply to it a pulse of an oscil- 
lating electric field with energy Wpinp=0.65 again. The 
time evolution of the ionicity p{t) is plotted in Fig. 5(a). 
For comparison, the dimerization (— l)-'wj(t) is shown in 
Fig. 5(b), and the displacement on the donor molecule 
—VT:i{t) in Fig. 5(c). Now, {—lyuj rapidly increases with 
p{t) and oscillates like (—1) times the cosine function. 
Thus, p{t) receives a positive feedback from (— 1)^-Uj from 
an early stage. It is also clear that the oscillations of uj™'' 



are temporally modulated by (— l)^Wj , as observed exper- 
imentally.^^^ In our calculation, the molecular vibrations 
and lattice phonons are indirectly coupled through elec- 
trons. 

In this calculation, fluctuations are not introduced 
into Uj, Uj, or wj™''. This leads to the fact that 

an ultrafast charge transfer between neighboring donor 
and acceptor molecules continues to oscillate without de- 
phasing. Although it appears as a very thick curve in 
Fig. 5(a), an ultrafast and large- amplitude oscillation is 
apparent if magnified on the time axis. In experiments 
performed at finite temperatures, such electronic motion 
is rapidly dephased. Indeed, if random numbers arc intro- 
duced according to the Boltzmann distribution at a finite 
temperature of 0.01 eV in Uj, v^"^\ iij, and uj™^ of the 
initial state, such an ultrafast and large-amplitude oscil- 
lation disappears but the slow oscillations caused by Uj 
and uj™^ survive. The resultant evolution of p{t) becomes 
closer to that observed in rcf. 24. These fiuctuations in- 
crease the ionicity on both sides of the phase boundary, 
which also becomes close to the experimentally observed 
behavior. As a consequence of the increased ionicity, the 
phase boundary is shifted to a larger A. We show an ex- 
ample of transient ionicity p(t) in Fig. 6 in such a case 
of dephased charge transfer in a neutral state near the 
shifted boundary. The ultrafast and large-amplitude os- 
cillation [the very thick curve in Fig. 5(a)] is indeed sup- 
pressed here. 

The above facts indicate the presence of strings of lat- 
tice dimerization before the photoexcitation of the neu- 
tral phase in the experiment performed at 90 K above 
the neutral-ionic transition temperature, i.e., 81 K.^''^ 
They are also consistent with the X-ray diffuse scatter- 
ing showing the characteristic length of 16 molecules for 
thermally induced ionic strings at 105 K.^s) At 90 K, 
the characteristic length would be larger and already be- 
yond the size for which we can treat the exact many- 
electron wave function. Such one-dimensional precursors 
locally break the inversion symmetry and make the quick 
growth and the — cosinc-type oscillation of the dimeriza- 
tion possible. Considering the fact that the photoinduced 
neutral-to-ionic dynamics decays much faster than the 
photoinduced ionic-to-neutral dynamics, ^''^ the growth 
of such one-dimensional ionic domains is saturated be- 
fore they form a three-dimensional ionic domain that is 
sufficiently metastable for longevity. 

The picture of the photoinduced neutral-to-ionic tran- 
sition obtained here is different from the previous one,^^^ 
which suggested ultrafast growth of an ionic domain 
without dimerization, based on the static and regular 
neutral initial state without couplings to lattice phonons 
or to molecular vibrations. In reality, the neutral phase 
in equilibrium above the phase transition temperature 
contains strings of lattice dimerization, already before 
photoexcitation. Photoexcitation rapidly grows and os- 
cillates the dimerization [Fig. 5(b)], which enhances the 
charge transfer between neighboring molecules [Fig. 5(a) 
compared with Fig. 4(a)]. In other words, precursors are 
transformed into ionic domains after photoexcitation. 

For clarity, the Fourier transform of p{t) during 50 < 
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t < 1550 [longer than the period shown in Figs. 4(a) 
and 5(a)], defined by | /^p e^'^*p{t)dt |, is plotted in 
Figs. 7(a) and 7(b). The presence of a low- frequency 
component at uj <0.02 in Fig. 7(a) is due to the slowly 
growing dimerization contributing to p{t) after the pe- 
riod shown in Fig. 4(a). When the dimerization is in- 
troduced in the initial state, the dimerization oscillates 
from the beginning [Fig. 5(a)] and contributes more to 
the low- frequency component in Fig. 7(b). With the ex- 
ception of such a low-frequency component. Figs. 7(a) 
and 7(b) are quite similar. All of the molecular vibra- 
tions evidently contribute to p{t): a peak at w=0.038 
due to VAi{t) with a;^Ai=0.040, a peak at w=0.050 due 
to vuit) with a;^D==0.055, and a peak at a;=0.11 due to 
VA2{t) with ci;^A2=0.12. All of them oscillate in the same 
phase as that in Figs. 4(b)-4(d). 

7. Conclusions 

Neutral-ionic transition and photoinduced dynamics 
have intensively been studied for the mixed-stack charge- 
transfer complex TTF-CA. The importance of the cou- 
pling to lattice phonons has been recognized from the 
dimerization and the consequent ferroelectric order in 
the ionic phase. The roles of the couplings to molecu- 
lar vibrations have not been paid much attention to, ex- 
cept for a few theoretical studies. ^^'"^^^ They are found to 
be important for the ionic phase to be a Mott insulator 
with large ionicity. The couplings to molecular vibrations 
stabilize the neutral phase, making the ionicity in the 
ionic (neutral) phase near the boundary large (small). 
The couplings to lattice phonons (molecular vibrations) 
reduce (enhance) the discontinuities in physical quanti- 
ties. Both couplings are necessary for the reproduction 
of the experimentally observed ionicity, optical conduc- 
tivity, and photoinduced charge dynamics. 

The photoinduced ionicity dynamics also shows a large 
contribution from molecular vibrations. The comparison 
with the experimentally observed, photoinduced charge 
dynamics indicates the presence of strings of lattice 
dimerization as local and precursory symmetry breaking 
in the neutral phase above the transition temperature. 
Only when a small but finite dimerization is introduced 
in the initial state can we reproduce its —cosine behavior 
and its large contribution to the photoinduced ionicity. 

Acknowledgment 

This work was supported by Grants-in-Aid for Scien- 
tific Research (C) (Grant Nos. 19540381 and 23540426), 
Scientific Research (B) (Grant No. 20340101) and Scien- 
tific Research (A) (Grant No. 23244062), and by "Grand 
Challenges in Next-Generation Integrated Nanoscience" 
from the Ministry of Education, Culture, Sports, Sci- 
ence and Technology of Japan, and the NINS program 
for cross-disciplinary study (NIFS10KEIN0160). 



1) Y. Tokura: J. Phys. Soc. Jpn. 75 (2006) 011001. 

2) K. Yonemitsu and K. Nasu: Phys. Rep. 465 (2008) 1. 

3) R. H. McKenzie: Science 278 (1997) 820. 

4) H. Seo: J. Phys. Soc. Jpn. 69 (2000) 805. 



10 
11 

12 

is: 

14 

is: 

16 
17 

18 

19 
20 

21 

22 

23 

24 

25 

26 
27 
28 
29 

30 

32 

33 
34 

35 

36 

37 
38 
39 
40 



Y. Tanaka and K. Yonemitsu: J. Phys. Soc. Jpn. 77 (2008) 
034708. 

S. Miyashita and K. Yonemitsu; J. Phys. Soc. Jpn. 77 (2008) 
094712. 

Y. Tanaka and K. Yonemitsu: J. Phys. Soc. Jpn. 79 (2010) 
024712. 

S. Miyashita, Y. Tanaka, S. Iwai, and K. Yonemitsu: J. Phys. 
Soc. Jpn. 79 (2010) 034708. 

Y. Kawakami, T. Fukatsu, Y. Sakurai, H. Unno, H. Itoh, 
S. Iwai, T. Sasaki, K. Yamamoto, K. Yakushi, and K. Yone- 
mitsu: Phys. Rev. Lett. 105 (2010) 246402. 
J. B. Torrance, J. E. Vazquez, J. J. Maycrlc, and V. Y. Lee: 
Phys. Rev. Lett. 46 (1981) 253. 

J. B. Torrance, A. Girlando, J. J. Maycrle, J. I. Crowley, V. Y. 
Lee, P. Batail, and S. J. LaPlaca: Phys. Rev. Lett. 47 (1981) 
1747. 

S. Koshihara, Y. Tokura, T. Mitani, G. Saito, and T. Koda: 
Phys. Rev. B 42 (1990) 6853(R). 

M. Le Cointe, M.-H. Lemeo-Cailleau, H. Caillcau, B. Toudic, 
L. Toupet, G. Heger, F. Moussa, P. Schweiss, K. H. Kraft, and 
N. Karl: Phys. Rev. B 51 (1995) 3374. 

S. Koshihara, Y. Takahashi, H. Sakai, Y. Tokura, and T. Luty: 
J. Phys. Chem. B 103 (1999) 2592. 

T. Suzuki, T. Sakamaki, K. Tanimura, S. Koshihara, and 
Y. Tokura: Phys. Rev. B 60 (1999) 6191. 

S. Iwai, S. Tanaka, K. Fujinuma, H. Kishida, H. Okamoto, 

and Y. Tokura: Phys. Rev. Lett. 88 (2002) 057402. 

E. Collet, M.-H. Lemee-Cailleau, M. Buron-Le Cointe, H. Cail- 

leau, M. Wulff, T. Luty, S. Koshihara, M. Meyer, L. Toupet, 

P. Rabiller, and S. Tcchert: Science 300 (2003) 612. 

L. Guerin, E. Collet, M.-H. Lemee-Cailleau, M. Buron-Le 

Cointe, H. Cailleau, A. Plech, M. Wulff, S. Koshihara, and 

T. Luty: Chem. Phys. 299 (2004) 163. 

K. Tanimura: Phys. Rev. B 70 (2004) 144112. 

H. Okamoto, Y. Ishige, S. Tanaka, H. Kishida, S. Iwai, and 

Y. Tokura: Phys. Rev. B 70 (2004) 165202. 

H. Matsuzaki, H. Takamatsu, H. Kishida, and H. Okamoto: 

J. Phys. Soc. Jpn. 74 (2005) 2925. 

G. D'Avino, A. Girlando, A. Painelli, M.-H. Lemee-Cailleau, 
and Z. G. Soos: Phys. Rev. Lett. 99 (2007) 156407. 

L. Guerin, J. Hebert, M. Buron-Le Cointe, S. Adachi, 
S. Koshihara, H. Cailleau, and E. Collet: Phys. Rev. Lett. 
105 (2010) 246101. 

H. Uemura and H. Okamoto: Phys. Rev. Lett. 105 (2010) 
258302. 

N. Nagaosa and J. Takimoto: J. Phys. Soc. Jpn. 55 (1986) 
2745. 

N. Nagaosa: J. Phys. Soc. Jpn. 55 (1986) 2754. 

A. Painelli and A. Girlando: Phys. Rev. B 37 (1988) 5748. 

Y. Toyozawa: Solid State Commun. 84 (1992) 255. 

T. lizuka-Sakano and Y. Toyozawa: J. Phys. Soc. Jpn. 65 

(1996) 671. 

P. Huai, H. Zheng, and K. Nasu: J. Phys. Soc. Jpn. 69 (2000) 
1788. 

Y. Anusooya-Pati, Z. G. Soos, and A. Painelli: Phys. Rev. B 
63 (2001) 205118. 

T. Kawamoto, T. lizuka-Sakano, Y. Shimoi, and S. Abe: Phys. 
Rev. B 64 (2001) 205107. 

K. Yonemitsu: Phys. Rev. B 65 (2002) 205105. 

L. Del Freo, A. Painelli, and Z. G. Soos: Phys. Rev. Lett. 89 

(2002) 027402. 

N. Miyashita, M. Kuwabara, and K. Yonemitsu: J. Phys. Soc. 
Jpn. 72 (2003) 2282. 

J. Kishinc, T. Luty, and K. Yonemitsu: Phys. Rev. B 69 (2004) 
075115. 

K. Yonemitsu: Phys. Rev. B 73 (2006) 155120. 

K. Iwano: Phys. Rev. Lett. 97 (2006) 226404. 

Z. G. Soos and A. PaineUi: Phys. Rev. B 75 (2007) 155119. 

K. Yonemitsu and N. Maeshima: Phys. Rev. B 79 (2009) 

125118. 



J. Phys. Soc. Jpn. 



Full Paper 



Author Name 



0.35 




la=0. 167, A,j=0.20,/1=0.219(N) 
£Bpmp=0.65,7tDmp=10,F=1.4 — 



100 200 300 400 500 600 
~1 




0.04 



0.03 



0.02 



0.01 



100 200 300 400 500 600 

t 



100 200 300 400 500 600 
t 



Fig. 4. (Color online) Transient quantities during and af- 
ter charge-transfer photoexcitation of neutral phase using 



=0.65, Tp, 



=10, and F=1.4 in case of strong Holstein and 



Peierls couplings Aa=0.167 and A^=0.20: (a) ionicity p(t) and 
displacements (b) VAi{t) with bare energy i.i;^ai=0-040 for the 
acceptor molecule, (c) —VY){t) with bare energy lj^d=0.055 for 
the donor molecule, and (d) v\2(.t) with bare energy ij^A2=0.12 
for the acceptor molecule. 



8 J. Phys. Soc. Jpn. 



Full Paper 



Author Name 



0.8 
0.6 
0.4 
0.2 


-0.2 



Aa=0. 167, A/j=0.20.4=0.219(N) 
topmp=0.65,7|Dmp=10.f=1.4 
(-l)'«,(r=-20)=0.01 
J '. I I I 



200 400 600 800 1000 



1000 




1 

0.8 k 
0.6 
0.4 
0.2 


-0.2 




AK=0.167,A,j=0.20,il=0.30(N) 
a)pmp=0.65,rpmp=10,F=4.2 
((-l>'«/r=-20)>=0.01 







200 400 600 800 1000 



t 



Fig. 6. Transient ionicity p(t) after setting initial dimerization 
{—l)^Uj(t = —20) = 0.01, adding random numbers to phonon 
variables as explained in text, and charge-transfer photoexcita- 
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Fig. 5. (Color online) Transient quantities after setting initial 
dimerization {—l)^Uj{t = —20) = 0.01 and charge-transfer pho- 
toexcitation of neutral phase using cjpinp=0.65, Tpmp=10, and 
F=1A in same case as that in Fig. 4: (a) ionicity p{t), (b) dimer- 
ization {—l)^Uj{t), and (c) displacement — i^dC*) with bare en- 
ergy aj^D=0.055 for the donor molecule. 
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0.01 set artificially. 



(b) 



A„=0.167,A,3=0.20,zi=0.219(N) 
fuext=0.65,rext=io,£ext=i.4' 
(-l)'M,(r=-20)=0.01 




